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Abstract. In this article, we collect the recent results concern- 
ing the representations of integers as sums of an even number of 
squares that are inspired by conjectures of Kac and Wakimoto. 
We start with a sketch of Milne's proof of two of these conjectures. 
We also show an alternative route to deduce these two conjectures 
from Milne's determinant formulas for sums of 4s^, respectively 
4s(s + 1), triangular numbers. This approach is inspired by Za- 
gier's proof of the Kac- Wakimoto formulas via modular forms. We 
end the survey with recent conjectures of the first author and Chua. 



1. Introduction 

The problem of finding explicit formulas for the number of represen- 
tations of an integer n as a sum of s squares is an old one. The first 
formula of this kind is due to Legendre and GauB. If rs{n) denotes the 
number of representations of n as a sum of s squares, Legendre and 
Gaufi proved that 

(1.1) r2(n) =4(rfi(r2)-d3H), 

where dj{n) denotes the number of divisors of n of the form 4k + j. 
For example, if is a prime p of the form 4k + 1, then r2{p) = 8 since 
di{p) = 2 and d^lp) = 0. On the other hand, if n is a prime p of the 
form 4k + 3, then r2{p) = since di{p) = d^i^p) = 1. This, of course, 
leads to the well known result of Fermat, which states that a prime p 
is of the form + if and only if p is of the form 4k + 1. Fermat 's 
result led mathematicians to explore and characterize primes of the 
form x"^ + ny'^, n > 1. For more information on such characterizations, 
the reader is encouraged to consult the excellent book by D. A. Cox 
0. 
Let 

oo 
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It is clear that 

fc>0 

As a result, to obtain expressions for rs{n), it suffices to obtain expres- 
sions for (p'^iq). The ffist identity of this kind is due to Jacobi, namely, 



— 2k— 1 



fc=l 

Note that (jl.lj) is a direct consequence of ()1.2j) . Using the theory of 
elliptic functions, Jacobi also found formulas for r4^{n), rQ^n) and rsin), 
namely, 

(1.3) /(,)=i+85:^-^ 



9) 

A:=l ^ fc=l ^ 



fc=l 



2k-l 



and 



(1.5) ^'(?)=i+i6f:Y^J^,. 

fc=l ^ 

From (|1.3p . we find that 



r4(r;,) = 8 



d^O (mod 4) 

and this immediately implies that every positive integer is a sum of 
four squares, a famous result of Lagrange. 
We call series of the type 



A 



q- 



q'' 

k>l ^ 



generalized Lambert series. Note that for even s < 8, we are able to 
express (p'^{q) in terms of generalized Lambert series. This does not 
seem possible when s = 10. In fact, Liouville showed that 



32 



k=l ^ k=l ^ 
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where 



oo 



fc(fc+l)/2 



fc=0 

Indeed for any even s > 10, </3*(q') is a sum of generalized Lambert 
series and a "cusp form" . 

Recently, new formulas for rs{n) were discovered. One common fea- 
ture of these formulas is the absence of "cusp forms''.^ The new formu- 
las involve only generalized Lambert series. The purpose of this article 
is to describe these recent discoveries. 

Before we proceed with our discussion, we make the following obser- 
vation: It is known that [21 p. 43, Entry 27 (ii)] 

(1.6) 4e"^/(2-)^2/ -2Wn = I^s.^^irN 

i 

Suppose we have a relation 

4.^./2^2s^g2) ^ F{L,{q''), L2{q^), . . . , LUq')), with q = e^'^ 

where each Lj is a generalized Lambert series or a product of general- 
ized Lambert series satisfying 

LJe-^^'/^) = (-y L* 



for some L*{—q) (which is also a generalized Lambert series or a prod- 
uct of generalized Lambert series), then we would have 

^'^{-q) = F{Ll{-q),...,L*^{-q)). 

Conversely if we have a formula for sums of squares, we will have a 
formula for sums of triangular numbers. We illustrate the above obser- 
vation by the following identities: 
Suppose for q = e'^*'^, we have 

7 3 2/c 

(1.7) 4V-V(e2--) = 16 5^^^ 

fc>0 ^ 

= ^ {E.ir) - E,{2t)) , 



where 



'1.8) E4(r) = 1 + 240"^ 



fc>l 



^Note that the coefficients of g" in Lambert series can be calculated once we 
know the factorization of n. In general, this is impossible for cusp forms. Hence, 
these new formulas are more "effective" if one wants to determine (n) ■ 
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The Eisenstein series -E4(t) satisfies the transformation formula PI p. 
24, Ex. 12] 



E,i--]= t^EAt). 



T ^ 

Replacing r by — l/r in ()1.7|) . the left hand side of ()1.7|) is T^^^{—q) 
by ()1.6|1 . with g = e"*'^. Now, by ()1.8|1 . we have 

- (i?4(-l/r) - i54(-2/r)) = r'- iE,{r) - -^E,{r/2) 



g,^ ' 1 — Q2fc+1 

fc>l ^ fc>l ^ 



Hence, we conclude that 

fc=i ^ 

Replacing q by — g, we obtain the formula for sums of 8 squares. For 
more details of the relations between formulas associated with squares 
and triangular numbers see [H] and PU]- 

Note that the identity for ip^{q) is much simpler than that for (p^{q). 
This is in fact a general phenomenon: the identity for ip'^'''{q) will be 
much simpler than that for (f'^'^{q) for any s G N. For the rest of this 
article, we will therefore only present identities associated with ip{q). 

2. The formulas of Kac and Wakimoto 
In 1994, V. G. Kac and M. Wakimoto [7] conjectured that 

(2.1) t4,2(n) = ^ ^2s-i ., ai ■ ■ ■ a^Ylia^i - a^jf 

■llj=l J- ai,--- ,aseN, ai odd i<j 
ri,--- ,rsSN, r; odd 

and 
(2.2) 

2 2^ -V -I— r 

tMs+i){n) = - ^ {ai---a,fY[{a^i -a^jY. 

■ 1 \.j=l3- ai,- ,asGN i<j 

'"ir" jJ^sSN, r'i odd 

These formulas follow from a conjectural affine denominator formula 
for simple Lie superalgebras of type Q{m). (For the definition of Q{m), 
see in].) 

Identities fj2.H) and ()2.2|) were first proved by S. C. Milne [H], us- 
ing results on continued fractions and elliptic functions. For example. 
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Milne showed using Schur functions, that (j2.Hl is a consequence of his 
determinant formula (5.107)] 

^2 4-s(s-i) 

(2.3) (g^^(g^)) = —2^^^det{C2(u+v-i)-i)i<u,v<s , 

n,=i 

where 

^2j-l - 2^ 1 _ 2(2r-l) ' ^ ^ 1- 

r=l 

We now briefly describe Milne's proof of ()2.3|) . 

Milne first showed that if sn(M) := sn(M, k), dn(M) := dn(M,k), and 
cn(M) := cn(M, k) are the classical Jacobi elliptic functions, then P 
(2.44), (2.68)] 

sn(u)cn(n) _ 1 ^ 22'"+2(-l)'"-i u^™"! 
^ ■ ' dn(u) ~ k2 ^ i^^^^ '"~'(2m-l)! 

^ ' m>l ^ ' 



where 



'^''™(2m-l)!^ 



(2.5) z = ip\q) and k^ = IGg-"^'^^'^ 



Milne then showed that 



°° sn(n)cnfu) 



(2.6) / ^^^^r^rr ^ e-/^dn 
dn(u) 



t2 



1 + (4 - 2k2)t2 + j^;;; 



oo -(2r2- l)(2n-2)2(2n-3)k^t^ 



1 + (2n- 1)2(4 -2k2)t 



2 



Here, KJ^g is the notation for continued fractions. 
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Using [HI Theorem 3.4] and ()2.fi|l . Milne deduced the Hankel determi- 
nant evaluation (4.9)] 



(2.7) i7«({c„}):=det 



C2 



C2 
C3 



Cn \ 

Cn+1 

C2n-l/ 



2n-l 



(1^2)n(n-l) -Q 



r=l 



Simplifying the left hand side of (j2.7p using the definition (j2.4j) of the 
Cj's and making use of the relations (3.66), (5.11)] 

and (EISl), we deduce (Q. 

We now describe a simplification of Milne's Schur function argument 
that allowed him to deduce ()2.H) from (j2.3|) . As a side result, we also 
obtain a new expression for t4s2(n) (see (12. 8p below). 

In a recent paper D. Zagier gave a direct proof of the above 
formulas of Kac and Wakimoto using the theory of modular forms. 
In that paper, he constructed a certain map sending the monomi- 
als X^^~^ ■ ■ ■ (here, the Xj's are indeterminates) to the product 
of Eisenstein series g'l^ - ■ ■ g'l^ (the quantities being, up to scaling, 
the quantities C2j-i in Milne's formula). It turns out that if we ap- 
ply a variant of that map, $s say, defined by sending the product 
■ ■ -Xf to the product C2fci-i ■ ■ ■C'2fc,-i, in Milne's formula, 
then we get the new formula 



(2.8) Us^{n) 



(_l)'^('^-l)/2 

ai,ri G N odd 



_ _3 _2s-l 



n 

l<i<j<s 



airiH |-asrs=2n+s^ 

This is seen as follows: By series expansion, we have 



2i-i 



r=l 



2fc(2r-l) 



k=0 



^^|2j-l^(2fc-l)(2r-l) 



r,fe>l 



m odd 
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Thus, we obtain 



^miH |-ms^2fci-l _ _ _ ^2fcs-l 



mi,...,ms odd 

ai\m-i, as|m. 



Returning to Milne's formula ()2.3p . this implies that 
det (C„„+„.„_0„„,,. = *, (det ,„,„^ 

= 4, fn^f-'detW-'>)_ 



-2^ 

1=1 l<*<j<s 



- a?) 



mi,...,ms odd i=l l<*<J<s 

ai|mi, as\ms 

where we have used the Vandermonde determinant evaluation to eval- 
uate the determinant in going from the second to the third line. Now 
a comparison of coefficients of g^n+s leads us to ()2.8|) . 

We now show that ()2.1|) follows from ()2.8p using an elementary com- 
binatorial argument. 

For each positive integer s, let 



and 



2^2 

j ■ 

j=l i<j 



p:{x„...,x,) = l[xr'l[{x^-x] 



-2^ 

j=i j<j 



For positive integers s and m, let 

i?,(m,P(Xi,...,X,)) = ^ P(ai,...,a,). 

o.i,ri G N odd 
Qir-iH \-asrs=m 

We want to show that 

Now, if Ss denotes the symmetric group of s elements and a & Sg, then 
P,(m, Pi(X,(i), . . . , = P,(m, Pi(Xi, . . . , X,)), 
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since, whenever (ai, . . . ,as) is an s-tuple of odd nonnegative integers 
for which there are ri, . . . such that aiVi + ■ ■ ■ + a^r^ = m, then 
{a„(i), . . . , ao-(s)) is an s-tuple with the same property. Hence, 

Rs (m, J2 P's{X.{ih • • - X,{s))) = s\ R,{m, P;(Xi, . . . , X,)). 

On the other hand, 

i?Ym, 5^Pi(X.(i),...,X^(,))') 

= i?Ym,det(Xf-2)i<,^.<,J^X, J] {Xl-X])\ 

by the Vandermonde determinant evaluation. Thus, we have shown 
how ()2.1|) follows from Milne's determinant formula ()2.3p by passing 
via ()2.8|) . thereby providing an alternative to Milne's (somewhat more 
involved) Schur function argument. 

3. A CONJECTURE FOR THE SUM OF 8s TRIANGULAR NUMBERS 

We now state Milne's determinant formula for 4s(s + 1) triangles: 
(3.1) (l6gV^^(g^))^^^^^) = (2^(-+^)) n(j!)-Met , 



where 

oo 



°° 2j+l 2r 



r=l 

This formula led to the first proof of ()2.2j) . Using arguments analogous 
to the ones given in the last section, one can deduce ()2.2j) from ()3.1|1 . 
When s = 2, this leads to the following beautiful formula: 



where 



n=l ^ 

Note the resemblance of this formula with the well-known formula 

00 



1728 

ri=l 
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where the i?j's are the classical Eisenstein series.^ Note that, as in- 
dicated in the introduction, one obtains Milne's new formula for 24 
squares, namely, if 



oo 



k'^q 



k=i ' 

and 



58(g) = 17 + 32^- 



fc=i 



-qy 



then P Theorem 1.6, (1.25)] 



^'\q) = l{S4{q)Ss{q)-8Sl{q)}. 

Comparing this with the "old" formula 

y,24^ ) ^ 1 — i?ii(g) + ^^qf\q) - ^^q'f\-q'), 
^ 691 691 691 ^ ^ ^ ' 

where 

oo 

/(-?)=n(i-^')' 

k=l 

we find that Milne's formula requires less terms. Moreover, if we know 
the factorization of n, then we can calculate r24,{n) explicitly from the 
new formula, since the terms are all Eisenstein series. The new formula 
for 24 squares and a recent paper of Z.-G. Liu [S] led the first author 
and Chua jH] to formulate the following conjecture for 8s triangular 
numbers: 

Conjecture. For any positive integer s > 1, we have 



m+n=2s 
m>n>2 



for some rational numbers Cr, 



For a fixed s one can verify the corresponding identity. For example, 
when s = 4 we are led to the following new identity: 

1 /95 91 \ 

(3.2) gV^(g2) = ^ [JT^o{q)Uq) - -T^{q) - Uq)Tu{q)j • 

^This was probably first observed by F. G. Garvan. 
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Note that the above identity does not follow from any of the formulas of 
Kac and Wakimoto or Milne, since 32 is not of the form or 4s(s + l). 

The first proof of this result proceeds by expressing the T2m's in 
terms of and z (see fl2.5|) for their definitions). The corresponding 
expressions are found by using the theory of modular forms, as well as 
the following recurrence satisfied by T^m^'- 

T2„+8(g) = T2(g)T2„+6(g) + ^^X] (2J + 2J ^2^+4 (g)T2„_ 2^+4 (g), 
where 



T2(g) = l + 24^- 



The above recurrence follows from the differential equation satisfied by 
the Jacobi elliptic function M := M[u) = sn^(u), namely. 



dM 
du 



4M(l-M)(l-k2M). 



We end this article with a sketch of a new proof of (j3.2p . It is known 
that the constant term of M{u) is T4 and that 

T4 = gV(g^). 

One can verify that 

(3.3) 3840M^ =(7M(^) + 68T4)(M(2) - 4T4) 

- 9M(=^)m(^) + M(2M(^) + 128T6), 

where M^*) is the ith derivative of M with respect to u. 

By comparing the constant term on both sides we are immediately 
led to ()3.2|) . The above identity is motivated by recent work of the 
first author and Liu j3] , where a proof of an identity similar to ()3.3|) is 
illustrated. 
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